Measurement of geo-neutrinos from 1353 days of Borexino 
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Abstract 

We present a measurement of the geo-neutrino signal obtained from 1353 days of data with the Borexino 
detector at Laboratori Nazionali del Gran Sasso in Italy. With a fiducial exposure of (3.69 ± 0.16) x 10^^ 
proton X year after all selection cuts and background subtraction, we detected (14.3 ± 4.4) geo-neutrino 
events assuming a fixed chondritic mass Th/U ratio of 3.9. This corresponds to a geo-neutrino signal Sgeo 
= (38.8 ± 12.0) TNU with just a 6 x 10"^ probability for a null geo-neutrino measurement. With U and 
Th left as free parameters in the fit, the relative signals are Sxh = (10.6 ± 12.7) TNU and S'u = (26.5 ± 
19.5) TNU. Borexino data alone are compatible with a mantle geo-neutrino signal of (15.4 ± 12.3) TNU, 
while a combined analysis with the KamLAND data allows to extract a mantle signal of (14.1 ± 8.1) TNU. 
Our measurement of 31.2tg ° reactor anti-neutrino events is in agreement with expectations in the presence 
of neutrino oscillations. 
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Geo-neutrinos (geo-Pe) are electron anti- 
neutrinos {i/e) produced mainly in /3 decays of '^''K 
and of several nuclides in the chains of long-lived 
radioactive isotopes ^ssy ^^^^ 232rp]_^^ which are 

naturally present in the Earth. By measuring the 
geo-neutrino flux from all these elements, it is in 
principle possible to deduce the amount of the 
radiogenic heat produced within the Earth, an 
information facing large uncertainty and being of 
crucial importance for geophysical and geochemical 
models. The first experimental investigation of 
geo-neutrinos from ^38^ g^j^j '^^'^Th was performed 
by the KamLAND collaboration [1] [2], followed 
by their observation with a high statistical sig- 
nificance of 99.997% C.L. by Borexino 3, and 
KamLANEj^ 4 . Both these experiments are using 
large-volume liquid-scintillator detectors placed 
in underground laboratories shielded against 
cosmic muons. Due to either low statistics and/or 
systematic errors, these measurements do not 
have the power to discriminate among several 
geological models. Analysis combining the results 
from different sites have higher prediction power, 
as it was shown in ^ and [6]. Therefore, new 
measurements of the geo-neutrino flux are highly 
awaited by this newly-born inter-disciplinary 
community. Several projects entering operation 
such as SNO+ [7 or under the design phase as 
LENA jSj or Hanohano [9j, have geo-neutrinos 
among their scientific aims. In this work we present 
a new Borexino measurement of the geo-neutrino 
signal with 2.4 times higher exposure with respect 
to [3]. For the first time, Borexino attempts a 
measurement of the individual geo-neutrino signals 
from the ^ss^j ^^^^ '^^'^Th chains. We provide a 
detailed comparison of our measurement with 
the predictions of several geological models. In 
a combined analysis of the Borexino and Kam- 
LAND 4J data we provide an estimate of the 
mantle geo-neutrino signal. 

Borexino is an unsegmented liquid scintillator de- 
tector built for the spectral measurement of low- 
energy solar neutrinos installed in the underground 
hall C of the Laboratori Nazionali del Gran Sasso 
(LNGS) in Italy. Several calibration campaigns 
with radioactive sources [10] allowed us to decrease 
the systematic errors of our measurements and to 



^An update of geoneutrino analysis from KamLAND col- 
laboration [5] was released after the submission of this paper 
and thus it was not considered in the combined analysis pre- 
sented below. 



optimize the values of several input parameters of 
the Monte-Carlo (MC) simulation. The 278 tons of 
ultra-pure liquid scintillator (pseudocumene (PC) 
doped with 1.5 g/1 of diphenyloxazole) are confined 
within a thin spherical nylon vessel with a radius of 
4.25 m. The detector core is shielded from external 
radiation by 890 tons of buffer liquid, a solution of 
PC and 3-5 g/1 of the light quencher dimethylph- 
thalate. The buffer is divided in two volumes by the 
second nylon vessel with a 5.75 m radius, prevent- 
ing inward radon diffusion. All this is contained 
in a 13.7 m diameter stainless steel sphere (SSS) 
on which are mounted 2212 8" PMTs detecting the 
scintillation light, the so-called Inner Detector. An 
external domed water tank of 9 m radius and 16.9 m 
height, filled with ultra-high purity water, serves as 
a passive shield against neutrons and gamma rays 
as well as an active muon veto. The Cherenkov 
light radiated by muons passing through the water 
is measured by 208 8" external PMTs also mounted 
on the SSS and define the so called Outer Detector. 
A detailed description of the Borexino detector can 
be found in milig. 

In liquid scintillator detectors, Ve are detected via 
the inverse neutron /3 decay, 

Pe + P — >■ e+ -I- n, (1) 

with a threshold of 1.806 MeV, above which lies 
only a small fraction of Pe from the 238]j (5,3%) and 
^•^^Th (3.8%) series. Geo-neutrinos emitted in "'"K 
decay cannot be detected by this technique. The 
positron created in this reaction promptly comes to 
rest and annihilates. All deposited energy is de- 
tected in a single prompt event, with a visible en- 
ergy of -Bprompt = Ei,^ - 0.784 MeV. The emitted 
free neutron is typically captured on protons, re- 
sulting in the emission of a 2.22 MeV de-excitation 
7 ray, providing a delayed coincidence event. The 
mean neutron capture time in Borexino was mea- 
sured with an AmBe neutron source to be r = 
(254.5 ± 1.8) /Lts [13]. The characteristic time and 
spatial coincidence of prompt and delayed events 
offers a clean signature of detection, further sup- 
pressing possible background sources. 

In this paper we report the analysis of data col- 
lected between December 2007 and August 2012, 
corresponding to 1352.60 days of live time. The 
fiducial exposure after all cuts is (613 ± 26) ton 
X year or (3.69 ± 0.16) x 10^^ proton x year. 

The Pe's from nuclear power plants are the main 
anti-neutrino background to the geo-neutrino mea- 
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surement. Since there are no nuclear plants close- 
by, the LNGS site is well suited for geo-neutrino 
detection. The expected number of events from re- 
actor i^e's, iVrcact, IS givcu by: 

R M 

Nreact = X! X! ^y^^™^ 
r— 1 m— 1 ^ 

(2) 

where the index r runs over the number R of re- 
actors considered, the index m runs over the total 
number of months M for the present data set, rjm is 
the exposure (in proton x yr) in the m*'^ month in- 
cluding detector efficiency, Lj. is the distance of the 
detector from reactor r, is the effective thermal 
power of reactor r in month to, the index i stands 
for the i-th spectral component in the set (^^^U, 
238u^ 239p^^ 24ip^)^ jg ^j^g power fraction 

of the component i, Ei is the average energy re- 
leased per fission of the component i, (j){Ei^) is the 
anti-neutrino spectrum per fission of the i^^ com- 
ponent, a{Ep) is the inverse /3 decay cross section 
taken from |14j . and Pee is the survival probabil- 
ity [6] of the reactor anti-neutrinos of energy Ep 
traveling the baseline Lr, for mixing parameters 9 
— {Sm'^, sin^6'i2, sin^^ia). 

In Eq. ([2| we consider the R — 446 nuclear cores 
all over the world, operating in the period of in- 
terest. The mean weighted distance of these re- 
actors from the LNGS site is about 1200 km, be- 
ing the weight Wrm = Prm/L'^- The effective ther- 
mal power, Prm, was calculated as a product of 
the nominal thermal power and the monthly load 
factor provided for each nuclear core by the Inter- 
national Atomic Energy Agency (IAEA) [T3]. For 
each core the distance, L^, has been calculated 
taking into account the geographic coordinates of 
the center of the Borexino detector (42.4540° lati- 
tude and 13.5755° longitude), obtained during the 
geodesy campaign for a measurement of CNGS 
muon-neutrino speed [IS]. The 0i(i?p_j) energy 
spectra are taken from [17j , differing from the spec- 
tra [18] used in [^ by about +3.5% in the normal- 
ization. The shapes are comparable in the energy 
window of our anti-neutrino candidates. Note that 
the 3.5% difference in the normalization is conser- 
vatively considered as a systematic error. For the 
power fractions, fi, we adopt the same assumptions 
as in our previous study 3J. Furthermore, in this 



analysis we precise fi for the 46 cores using heavy 
water moderator |19j : 

235u . 238^ • 239py . 241p^ _ 

(3) 

0.542 : 0.411 : 0.022 : 0.0243 ^ ' 

Since only two such cores are in Europe (in Ro- 
mania) this improvement in the calculation has an 
effect less than 0.1%. 

We adopt neutrino oscillations parameters as 
derived in |21j for normal hierarchy: Sm? — 
(7.54 tV^yiQ-^ eV2; sin^^i^ = (3.07 i°:l«)-10-i; 
sin^ Oxy. = (2.41 ± 0.25) -10-2. The three flavor sce- 
nario implies a 4.6% decrease in the predicted signal 
with respect the two neutrino case (as it was used 
in [3]), while the spectral shape does not signifi- 
cantly change. 

As in [3], we also include a -fO.6% contribu- 
tion from matter effects (oscillation parameters as 
above), and the -1-1.0% contribution of long-lived 
fission products in the spent fuel [50]. The contri- 
butions to the estimated systematic error are sum- 
marized in Tabled! 

Finally, the number of expected reactor De can- 
didates is Nreact — (33.3 ± 2.4) events for the ex- 
posure of (613 ± 26) ton x yr after cuts (for their 
efficiency see below). We note that in the absence 
of oscillation, the number of expected events would 
be 60.4 ±4.1. 

Table 1: Systematic uncertainties on the expected reactor- 
Pe signal which are added in quadrature. See Eq. |2| and 
accompanying text for details. 



Source 


Uncertainty 




[%] 


HE,) 


3.5 


Fuel composition 


3.2 


6*12 


2.3 


p 


2.0 


Long-lived isotopes 


1.0 




0.6 




0.5 




0.4 




0.4 


(5to^ 


0.03 


Total 


5.8 



The Borexino calibration campaigns [10| included 
several 7, /3, and a sources placed through the 
scintillator volume on and off-axis. The AmBe 
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source, producing ~10 neutrons/s with energies up 
to 10 MeV, was deployed in twenty-live different po- 
sitions allowing the study ol the detector response 
to captured neutrons and to protons recoiling off 
neutrons. The calibration data were essential for 
testing the accuracy of the Geant4-based Borex- 
ino MC simulation. The energy spectra of geo- 
neutrinos from ^ssy g^j^^j ^-^^Th, based on the theo- 
retical energy spectra of decays and the calcu- 
lated energy spectrum of reactor i?e's (see above), 
were used as input to the MC code in order to sim- 
ulate the detector response to Pe interactions. The 
MC output functions expressed in the total light 
yield, Q, (in units of photoelectrons, p.e., collected 
by the PMTs where 1 MeV corresponds to about 
500 p.e.) were then used as fit functions in the 
final analysis. In this way, the non-linearities of 
the detector response function important at higher 
energies and in the increased fiducial volume with 
respect to solar neutrino analysis, are automatically 
taken into account. 

The following cuts are used to select D^s can- 
didates: 1) Qprompt > 408 p.e. and 860 p.e. < 

Odelayod < 1300 p.C., whcrC Qprompt and Qdclayod 

are the PMTs' light yields for the prompt (positron 
candidate) and delayed (neutron candidate) events; 
2) reconstructed distance AR < 1 m; and 3) time 
interval 20 /iS < At < 1280 /iS between the prompt 
and the delayed event. In liquid scintillators, a 
pulse-shape analysis can be used to discriminate 
highly ionizing particles (a, proton) from particles 
with lower specific ionization (/3~, 7). The so- 
called Gatti parameter G [33] has been used to im- 
prove background rejection. For the delayed candi- 
date a very slight cut requiring Gddayed < 0.015 is 
applied. The total detection efficiency with these 
cuts was determined by MC to be 0.84 ± 0.01. 

A minimal distance of 25 cm from the inner vessel 
containing the scintillator is required for the posi- 
tion of the prompt candidate. Since this vessel is 
not perfectly spherical and does change in time, a 
dedicated algorithm was developed to calculate the 
vessel shape based on the position reconstruction 
of the events from the vessel's radioactive contami- 
nants. Since the vessel contamination is low, the 
vessel shape can be calculated only on a weekly 
basis. The precision of this method is 1.6%. It 
was calibrated by comparing the vessel shapes with 
those obtained by a dedicated LED calibration sys- 
tem |10j . The systematic error on the position re- 
construction of De candidates is 3.8% [3]. The total 
exposure of (613 ± 26) ton x year is calculated as 



a sum of weekly exposures which consider the cor- 
responding weekly live time and the vessel shape as 
well as the (0.84 ± 0.01) efficiency of the selection 
cuts described above. The 4.2% error on the ex- 
posure is a sum in quadrature of the errors on the 
vessel shape (1.6%), on the position reconstruction 
of the candidates (3.8%), and on the cuts efficiency 
(!%)■ 

Backgrounds faking anti-neutrino interactions 
can arise from cosmic muons and muon-induced un- 
stable nuclides, from intrinsic contaminations of the 
scintillator and of the surrounding materials, and 
from the accidental coincidences of non-correlated 
events. A complete list of all expected backgrounds 
is reported in Table |2] 

The levels of cosmogenic backgrounds {/3 + neu- 
tron decays of ^Li and ^He, fast neutrons, untagged 
muons) and of the background due to spontaneous 
fission in PMTs, did not change with respect to 
our previous paper 3 . We underline, that in order 
to suppress cosmogenic background we still apply a 
2 s veto after a muon passes through the scintillator 
(mostly for ^Li and ^He) and 2 ms veto after muons 
pass through only the water tank (mostly for fast 
neutrons). These vetos induce about an 11% loss 
of live time. In addition when possible, the pulse 
shape of the candidates was checked by an indepen- 
dent 400 MHz digitizer acquisition system in order 
to further suppress undetected muon background. 

The Borexino scintillator radioactivity has 
changed in time mostly because of the six purifica- 
tion campaigns performed in 2010 and 2011. Dur- 
ing periods of no operations, ^^"^Po, the main con- 
taminant important for i^e's studies, is observed to 
decay exponentially with a r = 199.6 days. The 
mean ^^"Po activity during the period used for this 
work is 15.8 counts/day/ton. Backgrounds from ac- 
cidental coincidences and from (a, n) interactions 
were evaluated according to the same methods as 
described in ^3 . 

During the purification campaigns some radon 
did enter the detector. The -^^^Rn has r = 5.52 days 
and within several days the correlated backgrounds 
disappear, leaving in the detector the correspond- 
ing amount of ^^'^Pb. These transition periods are 
not used for solar-u studies, but, with special care 
can be used for studies. The ^^''Bi(/3) - ^^'^Po(a) 
delayed coincidence has a time constant very close 
to the neutron capture time in PC. The a parti- 
cles emitted by the ^^^Po usually show a visible en- 
ergy well below the neutron capture energy window. 
However, in 1.04 x 10^** or in 6 x 10^^ of cases, the 
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^^''Po decays to excited states of ^^°Pb and the a is 
accompanied by the emission of prompt gammas of 
799.7 keV and of 1097.7 keV, respectively. In hq- 
uid scintillators, the 7 of the same energy produces 
more light with respect to an a particle [55] . There- 
fore, for these {a + 7) decay branches the observed 
light yield is higher with respect to pure a decays 
and is very close to the neutron capture energy win- 
dow. We have observed such candidates restricted 
to the purification periods, having the correspond- 
ing increased Qdciayod and positive (a-like) Gatti 
parameter. In order to suppress this background to 
negligible levels during the purification periods, we 
have increased (with respect to ^) the lower limit 
on Qdeiayed to 860 p.c. and applied a slight Gatti 
cut on the delayed candidate as described above. 



Table 2: Summary of the background faking anti— neutrino 
interactions and expressed in number of events expected 
among the 46 golden anti-neutrino candidates. The upper 
limits are given for 90% C.L. 



Events 



Background source 

^Li-^He 0.25±0.18 
Fast n's (^'s in WT) <0.07 
Fast n's (^'s in rock) <0.28 

Untagged muons 0.080±0.007 

Accidental coincidences 0.206±0.004 

Time corr. background 0.005±0.012 
(7,n) <0.04 
Spontaneous fission in PMTs 0.022±0.002 

(a,n) in scintillator 0.13±0.01 
(a,n) in the buffer <0.43 



Total 



0.70 ± 0.18 



We have identified 46 golden anti-neutrino can- 
didates passing all the selection criteria described 
above, having uniform spatial and time distribu- 
tions. All prompt events of these golden candidates 
have a negative G parameter, confirming that they 
are not due to a's or fast protons. The total number 
of the expected background is (0.70 ± 0.18) events 
(see Table [2]) . The achieved signal-to-background 
ratio of ~65 is high due to the extreme radio-purity 
of Borexino scintillator and high efficiency of the 
detector shielding. 

In the energy region Qprompt > 1300 p.e., above 
the end-point of the geo-neutrino spectrum, we ob- 
serve 21 candidates, while the expected background 
as in Table [2] is (0.24 ± 0.13) events. In this en- 
ergy window, we expect (39.9 ± 2.7) and (22.0 ± 



1.6) reactor-f/g events without and with oscillations, 
respectively. The expected survival probability is 
therefore (55.1 ± 5.5)%, a value almost constant for 
distances Lr > 300 km. We recall that for Borex- 
ino the closest reactor is at 416 km and the mean 
weighted distance is 1200 km. We conclude that 
our measurement of reactor D^'s in terms of num- 
ber of events is statistically in agreement with the 
expected signal in the presence of neutrino oscilla- 
tions. The ratio of the measured number of events 
due to reactor D^s with respect to the expected 
non-oscillated number of events is (52.0 ± 12.0)%. 

We have performed an unbinned maximal likeli- 
hood fit of the light yield spectrum of our prompt 
candidates. The weights of the geo-neutrino 
(Th/U mass ratio fixed to the the chondritic value 
of 3.9 [28j) and the reactor anti-neutrino spec- 
tral components were left as free fit parameters. 
The main background components were restricted 
within ±1(7 around the expected value as in Table[2j 
For the accidental background we have used the 
measured spectral shape, while for the (a, n) back- 
ground we have used a MC spectrum. For the ^Li 
and ^He background we have used a MC spectrum 
as well which is in agreement with the measured 
spectrum of 148 events satisfying our selection cuts 
as observed within a 2 s time interval after muons 
passing the scintillator. 

Our best fit values are iVgeo = (14.3 ± 4.4) events 



and A^react 

nals S„p„ 



84.5 



+ 19.3 
-16.9 



31.2j;g " events, corresponding to sig- 
(38.8 ± 12.0) TNt|^ and S-^oact = 



TNU. The measured geo-neutrino signal 



corresponds to overall I'e fluxes from U and Th de- 
cay chains of 0(U) = (2.4 ± 0.7) x 10*^ cm-^ 
and (/)(Th) = (2.0 ± 0.6) x 10^ cm'^ s-\ consid- 
ering the cross section of the detection interaction 
(Eq. [l]) from [H]. From the ln£ profile, the nuU 
geo-neutrino measurement has a probability of 6 x 
10~^. The data and the best fit are shown in Fig.[l| 
while Fig. H shows the 68.27, 95.45, and 99.73% 
C.L. contours for the geo-neutrino and the reac- 
tor anti-neutrino signals in comparison to expecta- 
tions. The signal from the reactors is in full agree- 
ment with the expectations of (33.3 ± 2.4) events 
in the presence of neutrino oscillations. 

A contribution of the local crust (LOG) to the to- 
tal geo-neutrino signal, based on the local 3D geol- 
ogy around the LNGS laboratory, was carefully esti- 
mated in 32J as S'gco(LOG) = (9.7 ± 1.3) TNU. The 



^1 TNU = 1 Terrestrial Neutrino Unit 
/ lO-^^ protons 



1 event / year 
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Light yield of prompt event [p.e.] 

Figure 1: Qprompt light yield spectrum of the 46 prompt 
golden anti-neutrino candidates and the best fit. The yel- 
low area isolates the contribution of the geo-Pe in the total 
signal. Dashed red line/orange area: reactor-Pe signal from 
the fit. Dashed blue line: geo-Pe signal resulting from the 
fit. The contribution of background from Tab. [2] is almost 
negligible and is shown by the small red filled area in the 
lower left part. The conversion from p.e. to energy is ap- 
proximately 500 p.e./MeV. 

contribution from the Rest Of the Crust (ROC), 
based on the recent calculation by Huang et al. [33] . 
results in the geo~neutrino signal from the crust 
(LOC-fROC) of Sgco (Crust) = (23.4 ± 2.8) TNU. 
Subtracting the estimated crustal components from 
the Borexino geo-neutrino rate, we can infer the 
contribution of the mantle, S'gco(Mantlc) = (15.4 ± 
12.3) TNU. 

On the basis of cosmochemical arguments and 
geochemical evidences, the different Bulk Silicate 
Earth (BSE) models predict the chemical compo- 
sition of the Primitive Mantle of the Earth sub- 
sequent to the metallic core separation and prior 
to the crust-mantle differentiation. The predicted 
amount of U and Th in the mantle can be obtained 
by subtracting their relatively well known crustal 
abundances from the BSE estimates. The mantle 
geo-neutrino signal on the Earth surface depends 
not only on the absolute abundances of the radioac- 
tive elements but also on their distribution in the 
present mantle. For a fixed mass of U and Th, the 
extreme cases of 5gco (Mantle) are obtained by dis- 
tributing their abundances cither homogeneously in 
the mantle (so called high model) or in an enriched 
layer close to the core-mantle boundary (so called 
low model) |34l I35j . In this perspective our results 
are summarized in Fig.|3j which is obtained by com- 
bining the expected geo-neutrino signal from the 
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Figure 2: The 68.27, 95.45, and 99.73% C.L. contour plots 
for the geo-neutrino and the reactor anti— neutrino signal 
rates expressed in TNU units. The black point indicates 
the best fit values. The dashed vertical lines are the la 
expectation band for Srea- The horizontal dashed lines show 
the extremes of the expectations for different BSE models 
(see Fig. [s] and relative details in text). 

crust (LOC + ROC) with those from different BSE 
models reported in Table V of 6 . The current re- 
sult cannot discriminate among the different BSE 
models. 

We have performed a combined analysis of our re- 
sult with that of KamLAND [4j in order to extract 
the ^geo (Mantle). First, the corresponding LOC + 
ROC crustal contributions taken from [5] and [55] . 
respectively, were subtracted from the measured 
S'goo signal: Sgeo (Crust) = (23.4 ± 2.8) TNU for 
Borexino and 5'geo(Crust) = (25.0 ± 1.9) TNU for 
KamLAND. Then, a spherically symmetric mantle 
was assumed. The best fit value for the mantle sig- 
nal common for both sites is S'goo (Mantle) — (14.1 
± 8.1) TNU. 

The Earth releases radiogenic heat, -ffgco, to- 
gether with geo-neutrinos in a well fixed ratio, 
however the observed geo-neutrino signal depends 
both on the abundances of the individual radioac- 
tive elements and on their distribution inside the 
Earth. To extract the radiogenic heat power from 
a measured S'goo is therefore model dependent. We 
have calculated the expected Sgoo (U-I-Th) as a func- 
tion of the radiogenic heat produced by U and 
Th, iJgoo(U-HTh), for the Borexino and KamLAND 
sites (see Fig. [4]), and compared it to the Borex- 
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Figure 3: Geo— neutrino signal 5gco in Borexino (solid line) 
with itlcr uncertainty (dashed lines) compared with the pre- 
dicted values. The ilcr band of 5geo (LOG + ROG) crustal 
contribution [6] is summed with Sgeo (mantle) according to 
seven BSE models: a) Javoy et al. I31| , b) Lyubetskaya and 
Korenaga 1301 . c) McDonough and Sun I29| . d) AUegre et 
al. [27], e) Palme and ONeil |26], f) Anderson |25], g) Tur- 
cotte and Schubert [24|. Red (blue) segments correspond to 
high (low) models obtained with two extreme distribution of 
U and Th in the mantle as described in the text, based on [6]. 
On the z-axis we show the total uranium mass predicted by 
each BSE model in the primordial mantle. 



ino and KamLAND |3] results. The allowed re- 
gions between the red and blue lines in the plane 
<S'goo(U-l-Th) and i7goo(U-fTh) contain models con- 
sistent with geochemical and geophysical data. For 
each total mass of U and fixed Th/U ratio, the max- 
imal geo-neutrino signal (red line) can be obtained 
by maximizing the radiogenic material in the crust 
and allowing uniform distribution in the mantle. 
Similarly, the minimal signal (blue line) is obtained 
for the minimal radiogenic mass in the crust with 
the rest concentrated in a thin layer at the bot- 
tom of the mantle. The expected signal from the 
crust is taken from Table V of [6]. We have chosen 
as a reference the BSE model from 29J, predict- 
ing that the silicate Earth contains m(U) — (0.8 ± 
0.1) X 10^'' kg with mass ratios Th/U = 3.9 and 
K/U = 12000. The green regions are allowed by 
the BSE model The arrow "Min" indicates 

the contribution of the crust only. The arrow for 
the fully radiogenic model indicates 39.3 TW: it as- 
sumes that the total Earth surface heat flux of (47 
± 2) TW [3H] is completely due to radiogenic heat 
from U, Th, and K. Taking the relative proportions 
from the BSE of [2^, we get that in a fully radio- 
genic Earth, U, Th, and K produce 19.1, 20.2, and 
7.7 TW, respectively. 

We have performed another unbinned maximal 
likelihood fit of our 46 golden candidates in which 
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Figure 4: The signal ^u+Th from U and Th geo-neutrinos 
as a function of radiogenic heat production rate i?u+Th in 
Borexino (top) and KamLAND (bottom). Details in text. 



the individual contributions from the ^ssy g^j^^j 
^•^^Th chains were fitted individually (see Fig. [5|, 
with all other fit details as above. The best fit 
values are A^xh = (3.9 ± 4.7) events and Au = 
(9.8 ± 7.2) events, corresponding to S'tii — (10.6 ± 
12.7) TNU and S'u = (26.5 ± 19.5) TNU and 
fluxes (above OMeV) of (/)(Th) = (2.6 ± 3.1) x 10^ 
cm-2 s-i and </>(U) = (2.1 ± 1.5) x 10^ cm'^ s"!. 
The 68.27, 95.45, and 99.73% C.L. contour plots of 
5'Th versus are shown in Fig. [6j Although our 
data is compatible within Icr with only ^as^j signal 
(and ^Th = 0) or only 232Th signal (and S^u = 0), 
we note that the best fit of the Th/U ratio is in 
very good agreement with the chondritic value. 

A geo-reactor with thermal power <30 TW and 
235u . 238^ _ Q . q 23 composition was suggested 
by Herndon [35] . It is assumed to be confined in the 
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Figure 5: Qprompt light yield spectrum of the 46 prompt 
golden anti-neutrino candidates and the best fit with free U 
(blue) and Th (cyan) contributions. The yellow area isolates 
the total contribution of geo-i^eS. Dashed red line/orange 
area: reactor-Pe signal from the fit. The contribution of 
background from Tab. [2] is almost negligible and is shown by 
the small red filled area. The conversion from p.e. to energy 
is approximately 500 p.e./MeV. 
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Figure 6: The 68.27, 95.45, and 99.73% C.L. contour plots 
of the Sxh £^rid Su signal rates expressed in TNU units. The 
black point indicates the best fit values. The dashed blue 
line represents the chondritic Th and U ratio. 

central part of the Earth's core within the radius of 
about 4 km [37] . We have produced MC spectra of 
the expected geo-reactor anti-neutrino. In a simi- 
lar unbinned maximal likelihood fit of our 46 golden 
anti-neutrino candidates we have added another fit 



component, iVgco-ioact, while constraining iVi-oact to 
the expected value of (33.3 ± 2.4) events. All other 
fit details were as above, including fixed chondritic 
mass Th/U ratio. We set the upper limit on the 
geo-reactor power 4.5 TW at 95% C.L. 
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